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PREFACE 

The purpose of this study was to fill a need of the Space 
Operations Directorate of the North American Aerospace Defense Command 
(NORAD). One of the missions of the Space Operations Directorate is 
to maintain positional data on all man made objects in space. Ocassionally 
some of these objects cannot be located. To find a lost object, it 
is neccessary to estimate where the object could be from historical 
positional data on the satellite. This is simple to do over a short 
time but the error increases rapidly over longer times because of changes 
in the satellite's period due to orbital perturbations. For near earth 
satellites one of the most important perturbations is atmospheric drag 


which is influenced by solar activity. A method to predict this change 


in period (or mean motion) from initial data was needed. Required 


characteristics for accomplishing this are simplicity and computational 
efficiency since many objects must be updated at any time. 

The method developed in this thesis is simple and rapid. Tests 
were performed on actual data to verify the model. 

I would like to thank all those who helped in the production of this 
thesis. First of all, I would like to thank my advisor,LtC Charles Ebeling, 
for his assistance and patience on this project. I would also like to thank 
Lt Terry Sparks of NORAD for his help in collecting the data needed to 
carry out this thesis. I have special thanks for Mary Browning and 
Pam McCarthy of the library for their help in locating some very obscure 


references. 
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‘ n Mean motion 
: n First time derivative of mean motion 
n Second time derivative of mean motion 
; N/2 Differance betx2en Brouwer and Kozai mean motion 
P Period of the orbit 
? P Change in period per day 
t Time 
r at Change in time 
SF Solar flux F10.7 value 
s Change in SF at time t 


Change in SF at time t-dt 
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Drag force 
Drag coefficient 


“4 


P Atmospheric density 
; A Satellite's cross sectional area 
v Relative satellite velocity 

z Acceleration due to drag 

r Total satellite velocity 

mn Satellite mass 
‘ 
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\ ABSTRACT 
This investigation derived a simple model to determine the change 
in mean motion over time when the actual values are unknown. A method 
was developed to include effects of solar flux by calculating an average 
value of r over 30 days. The model requires a knowledge of the mean 


motion for about 30 days before the time of interest to calculate this 


average. 


The analysis was done using BMDP on a CDC Cyber 6000 computer 


using element set data from actual satellites. 


This model does not attempt absolute accuracy, but is intended 


to be a method to quickly approximate a new mean motion when real values 


le are not available. A limitation of this model is the amount of 


historical data and analyst judgement which are required. 
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STATISTICAL MODELS FOR PREDICTING THE CHANGE IN 
MEAN MOTION OF A SATELLITE OVER TIME 


INCLUDING THE EFFECTS OF SOLAR FLUY 


I. Introduction 


BACKGROUND 

Satellites in near earth orbit (those at an altitude of under 
1000 km) show a loss of orbital period known as decay. This decay is 
caused by the drag on a satellite due to the upper atmosphere. It would 
be a simple matter to account for this decay if the drag were constant. 
Drag, however, is not constant. It is a function of atmospheric density. 
This density is, in turn, a function of altitude and solar flux 
(changes in solar flux cause changes in the atmospheric density at all 
altitudes). It can be concluded, therefore, that the decay of a satellite 
is a function of its altitude and of the solar flux. 

There are several ways to calculate and predict this decay. One 
method requires an accurate model of the upper atmosphere. The model 
must include both solar flux and altitude dependence. One such atmospheric 
model is the Jacchia model. Using this model and the known parameters 
of the satellite, it is possible, through astrodynamics, to calculate a 
very accurate near term orbit for the satellite. This method is very 
time consuming and requires up-to-the-minute knowledge of the satellite's 
position. It also decreases in accuracy as the position is predicted 
further into the future (1). 


A second method is to calculate the rate of change in mean motion 
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‘ (the number of orbital revolutions per day, known as n) of a satellite. 

; « 

: & This is called NDOT. If NDOT is known, then the mean motion of the 
satellite and its position (or “element set") can be predicted for any 

p time. This prediction is known as “propagation of the element set." 
NDOT is found by a least squares fit of actual data to two different 
astrodynamic models (known as the SGP model and the GP4 model). This 
produces a fairly accurate model of the satellite's orbit in far less 
time than a complete model such as the Jacchia requires. It is also less 
accurate than the Jacchia model for the same time, but it loses accuracy 
less quickly than the Jacchia. Over longer times the new n produced 
from the simple model is more accurate than that from an exact model. That 
is, the model is valid for several days instead of hours. Both of the 
models are in use at the North American Aerospace Defense Command (NORAD). 

Ce While other models have been tried, none of them have outperformed 

these two. In tests done by NORAD, the models which provide some 
improvement in accuracy do so at the expense of greatly increased 


computer run times (2-3). 


RESEARCH OBJECTIVE 

A method that is accurate over long delta times and requires no 
knowledge of the satellite's current position is needed for predicting n. 
It should not require an increase in computer run times or size, above 
that of the current methods. It is not a replacement for any of the 
methods currently in use. It is to be considered as an additional 


method for an area where existing methods are weakest. 
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Research questions include: 
1) Is there a statistical relation between n and solar flux? 
2) What is the best statistical relation between n and n? 


3) What is the best statistical relation between n and n? 


4) Is there a statistical relation between n and other orbital elements? 


SCOPE 
There are several types of possible errors in satellites orbits 
that must be allowed for. These errors may be divided into three main 
types. They are altitude errors, plane errors, and in-track errors. 
Altitude errors are errors in the satellites orbital altitude. Plane 
errors are errors in the satellites orbital plane and can include errors 
in inclination and errors in right ascension of node. In-track errors 
are time bias errors between the time a satellite should be at some place 
and the time it actually arrives there. Altitude errors are normally 
measured in units of length. Plane errors are measured in units of 
angle. In-track errors are measured in terms of time. These errors 
may be combined for a satellite and given in terms of absolute distance 
error, but that is not normally done in routine cases. This thesis will 
concentrate on in-track errors only. 
This thesis will develop and test some methods of predicting mean 


motion of an orbit. Chapter II will begin this development by providing 


some background and basic theory of orbital decay. It will define the terms 


that will be used in this thesis, and will explore some of the current 
models and methods used for satellite decay. Chapter III will then 


develop the models that will be tested in this thesis. Chapter IV will 
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study the test results of these models and compare them to each other 
and a current model. Chapter V will then present the conclusions of this 
study. Extensive use will be made of tables and figures to illustrate 


the study that was done. 
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II. THEORY 


ORBITAL PARAMETERS 

A satellite's orbit is described by several parameters. There 
are several different sets of parameters used for this. These sets 
include the set of position and velocity vectors, the Keplerian element 
set, and the F& series (3). The parameters used in this thesis are 
those used by NORAD, which are almost the same as those in the Keplerian 
element set. There are some added parameters which will be included in 


the following discussion. 


The parameters used in the NORAD element set (elset) are: 


Epoch time: As used by NORAD in a general perturbations element 
set the epoch time of a satellite is taken as the time of passage 
through the equatorial plane on the last ascending pass of the satellite. 

Inclination: Inclination is defined as the angle made between the 
satellite's orbital plane and the equatorial plane in the direction of 
satellite movement at the ascending node. 

Ascending node: The ascending node is the point of passage of the 
satellite's orbit through the equatorial plane in a south to north 
direction. It is measured in degrees right ascension for the first 
point of Aries. 

Eccentricity: Eccentricity is a measure of the flatness of an 
orbit since the orbit is in the shape of an ellipise. 

Argument of perigee: The argument of perigee is the angular distance 
between the ascending node and the perigee of the orbit in the direction 


of satellite motion. 
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Mean anomaly: Mean anomaly is a measure of the angular distance 
between the argument of perigee and the satellite's position in the orbit, 
in the direction of satellite motion, at the epoch time. In the NORAD 
general perturbations element set the mean anomaly and argument of 
perigee must add to 360 degrees. 

Mean motion (n): Mean motion is a representation of a satellite’s 
orbital period. It is not, however, expressed in the time to complete 
one revolution, but in terms of revolutions per day. Its relation to 
period is given by n=1440/P, where P is in minutes. 

NDOT /2 (N/2)s NDOT/2 is a term given with the NORAD elset as the 
difference between a Brouwer and a Kozai mean motion (1), where Kozai 
used a fourth order general perturbations model and Brouwer used a 
simplified model of a satellite's orbit. It estimates orbital decay rate. 

BSTAR: BSTAR is a synthetic drag term for use in the equations of 
motion. It is a comparison of the orbit of a satellite to some standard 
reference satellite. It does not account for size or shape differences 
or solar flux. Its prime use is to indicate which satellites have been 
more affected by changes in solar flux. There is little other practical 
use for this tern. 

Figure 2.1 shows some typical NORAD 2-card elsets. 

Three other parameters are used in this thesis, though neither is 
in a normal elset. One of them is PDOT, which is the rate of change of 
orbital period per day. It is provided by the Naval Space Surveillance 
Center (NAVSPASUR). Another parameter used in this thesis is NDOT Gals 
For the purpose of this thesis n will be estimated as the difference 


between n at two different times, divided by the difference in time 


2.2 


such that: 


n(t.) = n(t,) - n(t,) (2.1) 


ty = ue 


The final parameter is n, which, in this thesis, is the difference 


in n at two times divided by the delta time so: 


n(t,) 2 n(t,) - n(t,) (2.2) 


17 ue 


The parameters N/2, n, n will be used in this thesis to estimate 


new values of n. 
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Figure 2.1 NORAD 2-card element set 
The Format of the 2-card elset is: 
Line 1: line number, satellite number, international designator, 
le epoch time, N/2, N/6, Bstar, elset number. 
Line 2: line number, satellite number, inclination, right ascension, 
eccentricity, argument of perigee, mean anomaly, mean motion. 
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ORBITAL MECHANICS 


A satellite in near earth orbit is affected by the upper reaches 
of the atmosphere. The satellite has a velocity of abcut 7 km/s relative 
to the upper atmosphere. As in the movement of any body through a 
medium, there will be a retarding force on the body opposite to the 


direction of motion. This retarding force is called drag and is given 


by: 


D=- Cpav- (2.3) 


where D = drag force 


C = drag coefficient 


p = atmospheric density 
A = satellite cross sectional area 
V = relative vehicle velocity (3-21) 


This will give an acceleration of the satellite such that: 


¥ = -_CApvr’ (2.4) 
2m 

where ? = acceleration due to drag 

m = mass of the satellite 

Z = total satellite velocity (4-423) 

These equations can be combined with the standard equations of 
motion and integrated to find the satellites orbit. The largest obstruction 
to this is finding the correct value(s) of p to use in the model. Both 


altitude and solar flux affect p and must be accounted for. 


2.5 


CURRENT MODELS 

a“ Satellite decay is covered by several models. These can be either 
estimating type models such as that developed in this thesis, or 
atmospheric density models which compute values of p for use in the 


equations of motion. Each type of model has advantages and disadvantages. 


re If accuracy is required and values of solar flux are known, then the 
atmospheric density models are best. If solar flux values are unknown 
or long time estimates are needed, estimating models are more suitable. 
The exact choice of model must be made in view of the product desired. 

Of the atmospheric density models, one of the most useful to date 
is the Jacchia atmospheric model. There are several versions of the 
Jacchia atmosphere model. One of them, the J65 model produces the 
density (p) values through a table “look up" method. It uses the values 
. ie of the solar flux to find density values ina table. It relates values 
of the solar flux at a wavelength of 10.7 cm to density given ina 


table compiled by Nicolet (2-B14). It should be noted that the 10.7 cm 
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flux (F10.7) does not heat the upper atnosphere. The upper atmosphere 
is heated by the extreme ultra violet (EUV) radiation which cannot be 


observed at the earth's surface, whereas the F10.7 flux can be measured 
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(2-B16). The J70 model uses both an average value of F10.7 taken at a 
time of t - 400 days and the changing daily value for a time lag of one 


, day such that: 


T = 383 + 3.3F + 1.8(F - F) (2.5) 


2.6 
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where F = average value of F10.7 flux at t-400 days 
ey F 


previous days value of F10.7 flux 


T = temperature (2-B17) 
The J70 model can then calculate the density through suitable manipu- 
lation of thermodynamic gas laws. 

Other models produce the p values in different ways. For example, 
the DENSEL model uses a power function in altitude and theoretical values 
of F10.7 (2-Bi4,B16) to find p. There is also the exponential atmos- 


phere: 


p= py exp-[e (h - h,)] (2.6) 


where Py = density at height hy 


a constant (2-B13) 


Q 
" 


ie The exponential atmosphere is similar +o the model used by Desmond 


King-Hele to determine satellite lifetimes, where King-Hele uses: 


p= p, exp [-(y - y)/t] (2.7) 


where B= density at Ye 


H scale height (5-182) 


and then: 
z = ae/H (2.8) 


where a = semi major axis 


e = orbital eccentricity (5-182) 


so that the satellite lifetime L may be written as: 


2.7 
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L=3en}1+7etSet_1/l+1le+ 3+ 3 o10.5/9) 
A 4al 6 16 2z\ 12 424 
¥ (5-182) (2.9) 


Any of the models that calculate density are accurate for the 
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time calculated, but are valid for only a short time before the orbit 
changes beyond prediction limits. This is due to uncertainty of future 


solar flux, and, therefore, atmospheric density values. Another problem 


ta. 


with this type of model is the calculation time required, on the order 


of 31 seconds to produce the densities (2-3), and minutes to hours to 


integrate the equations of motion for each satellite. This quickly 


grows beyond acceptable limits when dealing with many satellites. 


A method for predicting long term changes is in use by Naval Space 


Surveillance (NAVSPASUR). This method is based on calculating values of 


PDOT. PDOT is the change in period per day of a satellites orbit. It 


Le is not highly accurate, but offers two advantages: 


1) It gives a general indication of a satellites period over 


very long time spans. 


2) It is very quick, requiring less than a minute on a hand- 


held calculator. The equation used by this method is: 


P(t) = P(0) + Pat (2.10) 


This method is used by NAVSPASUR and NORAD to determine if a lost 


satellite may have decayed. Because of the inaccuracies, the satellite 


is not counted as decayed unless the period has been below 87.5 minutes 


(the standard cutoff for an orbit) for over a month, and there are no 


possible unknown satellites between its old position and possible 


current positions. In practice the use of this method is more an art 


than a science. It is only used by the most qualified analysts because 
a great deal of judgment on the part of the analyst is required. 

A third method, used mostly by NORAD, is based upon NDOT. NDOT 
is the difference between a Brouwer and a Kozai mean motion for the same 
time (1). Brouwer and Kozai are two methods of calculating mean motion 
used in the NORAD system. The method is very similar to the PDOT 
except that the change is calculated for mean motion instead of period. 


In application the equation used is: 
n(t) = n(0) + 2(N/2)at (2.11) 


Like the PDOT, the NDOT requires experience to use and produces 
similar results. 

The primary model developed in this thesis is based upon the NDOT 
model. It is employed in much the same manner as described above but 
includes an allowance for solar flux and error boundries on the accuracy 
of the calculated value of n. This provides a measure of the uncertainty 


in predicting a new value for n. Solar flux is addressed by calculating 


an average N/2 from several points out of a 30 day interval since solar 


flux varies over about a 27 day interval (6-189). 


MEASURES OF EFFECTIVENESS 

The accuracy limits are determined to allow a reasonable prediction 
at 100 days and still have a recoverable satellite. From experience 
this was set at about 0.1 min/rev which gives an error limit of 0.01595 
rev/day (for a 95 minute orbit). This is an error range of 1.58 revs at 


100 days. 
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III. METHODOLOGY 


DATA 

The data used in this thesis was provided by NORAD/HDS. The data 
consisted of 50 days of element sets on a group of 5 random satellites 
chosen from a set of 70 satellites with mean motions between 16 and 14 
rev/day. The satellites were chosen at random and were not given 
individual identifiers to eliminate the possibility of selecting data 
with some bias. The following data on each of the 5 satellites was 
collected: time, delta time, n(t), n(t-dt), n(t), n(t=0). 
In addition to these data, the values of F10.7 solar flux for the epoch 
date, the change in F10.7 flux from the previous day and the change in 
F10.7 flux on the previous day were added. This information was placed 
ina file on the Cyber computer. Appendix A contains a listing of this 
file. The data are in the format: 

t, dt, n(t), n(t-dt), n(t), n(O), S, F, SF. 

The values of n were calculated by equation 2.1 and placed in the 
file along with the data transferred from the NORAD elsets. 

The data in this file were then used in the program BMDP to fit the 


models to the data. 
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MODEL DEVELOPMENT 
A model for n may be developed in several ways. One of the models 
used at NORAD involves a Taylor series expansion of n to give: 


n(t) = n(0) + n(t - t,) + in(t - a +n% 
2 (3.1) 


where NORAD defined n as 2(N/2) and fi as 3(N/6) so that: 


n(t) = n(0) + 2(N/2) (t - t,) + 3(H/6) (t- t,)? +. . 
(3.2) 


and t - te is the time of interest minus the initial time. The second 
and higher order terms are of small magnitude compared to n for the 
normal values of (t - t,) encountered in element propagation. Therefore 
the second and higher order terms are set to zero. In normal usage the 


equation used for n is: 
n(t) = n(0) + 2(N/2) (t - t,) (3.3) 


That this formula is appropriate can be seen from figure 3.1 which 
is a plot of mean motion versus time. This figure shows that the value 
of n at any one time depends on the value of n at some previous time plus 
the sum of all the changes in n over the time difference. The sum of 
all changes can also be considered as an integral over all changes. 

Under the assumption that all changes over time greater than second order 
are very smz1l and can be considered zero, the second order change, n, 


can be considered a constant. That is: 


n=a (3.4) 
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Figure 3.1 Typical Mean Motion for days 84300 through 85084 
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% “a And since: 
es 
‘ ae an = an (3-5) 
¢ at~ at 
if 
; This gives: 
dn=a (3.6) 
: dat 
: Which may be written as: 
g dh = adt (3.7) 
: Equation 3.7 may be integrated: 
fan = afat (3.8) 
Le n-n, = a(t - t,) (3.9) 
Be 
LS 
“, a . 
; nein, a(t - t,) (3.10) 
: And since: 
° A = dn (3.11) 
of dt 
; dn =n, + a(t - t,) (3.12) 
dat 
5 dn = E + a(t - *_)|ae (3.13) 
ae 
a 3.4 
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Integrating equation 3.13 gives: 


ty: 
eed . 
S fan =ffn, + a(t - t Jat (3.14) 
psa See ee) tat eye (3.15) 
° ° fe) 3 ro) . 
‘ 2 
ee Ay(t - t,) + at - t.) (3.16) 
Substituting in equation 3.4 lead to an equation of the form: 
n(t) = n(o) + a(t ~ t,) + 18(t - t,)? (3.17) 
2 
This equation is basically the same as that given by a Taylor series, but 
is arrived at from physical considerations. 
From these come two of the models that will be studied in this ! 
le thesis. Equation 3.17 is the baseline nO model and is used in BMDP 
in modified form. Equation 3.3 is the standard NORAD model (herein 
called model n5). 
Model nO appears in this thesis in three basic forms. It appears 
in its original version as given by equation 3.17. There are also two 
modified versions which are used in BMDP. They include model ni which 
included regression coefficients in each term and has the intercept defined 
at Byn (0) such thats 
e Leda 2 
ni: n(t) = Byn(0) + Bun(t - t_) + 0.5B,n(t - t_) 
1 2 ie} 3 fe) (3.18) 
where the B's are the regression coefficients from BMDP. Next, under 
the assumption that the regression coefficients contain the physical 
Poe 
eae 
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constants and the solar flux dependence, and that long term average 
solar effects are contained in the intercept since over time flux 


approaches a constant gives: 


n2: n(t) = B, + Byn(0) + BpA(t - t,) + BAA(t - ee 


(3.19) 
Next, assuming there is no (t - tye dependence, that is, the 
satellites "forgets" earlier changes such that there is no second order 

time dependence, gives model n3: 
n33 n(t) = B+ Byn(0) + Bon(t - t,) + Bafi(t - t)) ea 
The final model comes from assuming a simple linear relation 


between the terms. Model n4 is: 


nh: n(t) = B, +B (t - t.) +B n+ B,n(0) + BF 


(3.21) 


Note that only model n4 contains solar flux explicity, where it is 
given as the change in the flux for time t-1. Solar flux changes are 
contained implicitly within the regression coefficients and the n and 
n terms in models ni, n2, and n3. The reason for this can be seen by 
comparison of figures 3.1 and 3.2, where figure 3.2 is a plot of solar 
flux for the same time frame as figure 3.1. It is obvious from the 
graph that the daily solar flux value has little long term effect on 
individual values of n. 

Likewise, a comparison of figure 3.2 and figure 3.3 (which is a 
graph of the daily change inn, i.e. n) shows a much higher correlation 
between n and flux. That this is indeed the case will be shown by the 


data analysis carried out in chapter four. 


ise we 


None of the models contains orbital elements other than those 


related to mean motion. Eccentricity is not included because it is 


near zero for the orbits being studied. While it is known that atmos- 
Y 
‘ pheric density is not constant around an orbit (density changes with 
: latitude and darkness), it is assumed in this thesis that effects re- 
, lating to inclination and other elements are very small for time spans 
: on the order of 100 days. Their inclusion in a model is an area for 
" additional study. 
. Below is given a complete list of the models that will be used 
in the analysis by BMDP, where t is now defined as t - t,. 
= : 42 

nd: n=n + nt + 0.5nt (3.17) 

: wo42 
: nis n= Byn, + Bont + 0.5B,nt (3.18) 
le 
z . wt 2 

n2: n= B+ Byn, + Bont + B,nt (3.19) 

nj: n= B+ Bun + Bont + B,nt (3.20) 
: ny; n= B+ Byt + Bon + Bn, + BF (3.21) 
, n5: n=n + 2(N/2)t (3.3) 
3.7 
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Figure 3.2 Solar Flux for days 84300 through 85084 
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Figure 3.3 Daily Change in Mean Motion for days 84300 through 85084 
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IV. RESULTS 


SOLAR FLUX 
In chapter three it was stated that solar flux was not a major 
explicit factor in predicting values of n. Some justification, on 
the basis of graphical analysis, was given. Further empirical justi- 
fication will be given here. 
In the BMDP analysis of model n4 it is shown that there is only 
a 0.0122 correlation between F and n (see Table 4.4), while the regres- 
sion coefficient is of the order of 10E-4. This shows that flux has 
little effect on n directly. A separate analysis shows that there 
is a -0.1021 correlation between F and n. 
Also, if flux is considered as a cumulative average and the daily 
. (a change from that average, it is seen that the average rapidly approaches 
a constant plus or minus the daily change. If this constant is included 
in the regression equations, it must be in the value of the intercept 
constant. The daily change from a constant will be in the daily values 
of the second and higher order terms of the n equation. That is, 
since for a nearly circular orbit there is no altitude dependence for 
drag, the rate of decay will contain a constant value for the average 
density at that altitude. This constant is the first order tern, n. 
The rate of change of n, in will be constant for a constant Change in 
solar flux. Since flux is not constant on a daily basis, there will be 
; a change in n, given by the third and higher order terms, which are 
nearly zero and are defined as zero such that n is a constant of small 


magnitude. Therefore the statement that the solar flux dependence is 
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implicit in n and the regression coefficients is supported. 


BMDP ANALYSIS 


The models included in the BMDP analysis were: 


ht + 0.5Bynt” 


ni: Byn +B 
° 


S 


. oe 2 
B, + Byn,, + Bont + B.nt 


nt + Bont 


B, + Byun, + By 3 


m4: B, + B,t + Bon + B,n,, + BUF 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


BMDP was run against the data file in Appendix A using each of 


these models independently. The output from BMDP produced the results 


shown in tables 4.1 through 4.4. 


2 NT 


ae TABLE 4.1 
b Re 
' BMDP RESULTS FOR MODEL N1 
Model nl: n= Byn, + Boat + 0. SBgit” (3.18) 
By aad ae EN cae ce ae a amet ee at 
Multiple R=1 R°=1 standard error=0.0118 
» B,=1 BA =1.06457 B,=0.01970 
: n=n +1 .06457nt+(0.0197/2) nt” 
Sum of Squares DF Mean Square Ratio P 
: Regression  22260.2993 3 7420.0998 53696658.325 0.0 
Residual 0.0130 94 0.0001 
Ce ee Se pe Swe te eg 
n, n n nt nt? 
Correlation--------------------------------------------- 
a with n 0.9962 0.8274 0.1138 0.7135 0.1118 
De. Re 
é 
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TABLE 4.2 
«,% 
oa BMDP RESULTS FOR MODEL N2 
: i panne 
Model n2: n= 8B + Byun + Bont + B,nt (3.19) 
Multiple R=0.9992 R“=0.9983 _tandard error=0.0115 
B=-0.20596 B,=1.01364 B,=0.9839u B,=0.00935 
n=-0,20596+1 .01384n, +0.98394nt+0. 0093 srt? 
Sum of Squares DF Mean Square Ratio P 
Regression 7.3171 3 2.4390 18347.738 0.000 
" Residual 0.0124 93 0.0001 
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TABLE 4.3 
BMDP RESULTS FOR MODEL N3 


One ww wm ew we ww oe ew we we ee www a ww ow ow ww oO nw we ww on we ew wm = ow = ee ew ee ew ww ee eee = 


Model n3: n= Bo + Bun + Bont + Bynt (3.20) 


Multiple R=0.9992 - R“=0.9983. ~— standard error=0.0114 
B =-0.18935 B,=1.01273 B,=0.99835 B,=0.43515 
n=-0,18935+1 .01273n +0 .99835nt+0.4351 5nit 


wwe we ee we ee ww em ew ew ee ww ow = on a a a wn we ww a ww we ww ow ow we eee ee ee wee = 


Sum of Squares DF Mean Square Ratio P 
Regression 7.3174 3 2.4391 18659.405 0.00 
Residual 0.0122 93 0.0001 

ny on n nt nt 
Correlation-------~----------------~------~------- 
with n 0.9962 0.8274 0.1138 0.7135 0.1089 
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TABLE 4.4 


BMDP RESULTS FOR MODEL N4& 


Model ny: on = B+ Byt + Bon + Ban, + EF (3.21) 


Multiple R=0.9981 R“=0.9962 standard error=0.0174 
B =-0.67573 B,=0.00090517 B,=17.45671 B,=1 . 04379 
B, =-0.009097511 


n=-0.67573+(9.0517E-4)t+17.45671 ntl .04379n +(-9.7511E-5)F 


www wm me we ewe ww www mw ewe ww wm wm we ww ww wm an we wm ww ww ww ww ww we we ww ww we eww wm em ew wwe eww we = 


Sum of Squares DF Mean Square Ratio P 
Regression 7.7315 4 1.9329 6366.816 0.000 
Residual 0.0294 97 0.0003 

t n n F 
Correlation------------------------ Beech Saas 
with n -0.0023 0.8160 0.9963 0.0122 


Pr 


wes 


= aT we . oy ee, ee 
wwe & fe af. ae A ort ay al Ma tee ee hae s 


POST BMDP ANALYSIS 

The equations produced by BMDP, along with nO and n5, were used 
to generate predicted n values for satellites 15363, 14476, 13043, and 
7840. The predictions were done for time spans of 10, 30, 50, and 100 
days (as close to these values as possible given the NORAD data). The 
calculated value of n was subtracted from the actual value of n at each 
time to give a value for the error inn such that a negative error 
indicates that the calculated value was greater than the actual value 
at time t. Tables 4.5 through 4.11 contain the sample data for these 
four satellites for each model. The values on n and n were calculated 


by equations 2.1 and 2.2, and are shown in Table 4.5. 


TABLE 4.5 


VALUES OF n, n, AND N/2 FOR EACH SATELLITE 


SATELLITE n n N/2 

15363 0.00061318 0.00007924 0.00033037 
14476 0.00002221 0.00000051 0.00001155 
13043 0.00049423 0.00001846 0.00027705 
7840 0.00002121 0.00000222 0.00001401 
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TABLE 4.6 


qn MEAN MOTION VALUES FOR MODEL NO 
SATELLITE : t ACTUAL n CALCULATED n an ERROR 
days rev/day rev/day rev/day 
(303, | Pico S206 ena ceeneon oO nee cco e co 
10 14.62055640 14.6260108 -0.00545440 
30 14.62850624 14..6699376 ~0.04143136 
50 14. 63611823 14.7455305 -0.10941227 
100 14.65988312 15.0730518 -0.41316868 
14176 “ibe Hipepeuogigy -\ceaeenacey ene cacee= 
10 14.47959223 14.4796510 -0.00005877 
31 14.48005842 14..4803389 -0.00028048 
52 14.48061500 14.4812537 -0.00063870 
102 14..48162641 14.4843450 -0.00271859 
13043 Wages, Meuegnediieg. ° eacucescae’ © get eae 
10 15.66745489 15.6687299 -0.00128501 
30 15.67900804 15.6860070 -0.00699896 
50 15.68914747 15.7106564 -0.02150893 
100 15.71437096 15.8045781 -0.09020714 
7840 0  14,56440052. ------=--- eee 
10 14. 56459216 14..5647237 -0.00013154 
Ke 29 14..56491839 14.5659499 -0.00103151 
52 14.56555106 14.5685076 -0.00295654 
103 14. 56649157 14.,5783721 -0.01188053 
4.8 
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TABLE +.7 


oS 
<> MEAN MOTION VALUES FOR MODEL N1 
SATELLITE t ACTUAL n CALCULATED 5. ERROF. 
days rev/day rev/day rev/day 
15363 0 Te61592219 eases CC ieee 
10 14.62055640 14.6225264 -0.00197000 
30 14.€2850624 14.6362026 -0.00769626 
50 1463611823 14.6505026 -0 05286367 
100 14,65988312 14.6889819 -0.02909878 
fee Pe ee Sie tena 
10 14.49959223 14.4796401 -0.00004787 
31 14.48005842 14.4801410 -0.00008258 
52 14.48061500 14.4806443 -0.00003130 
102 14.48162641 14.4818675 -0.00024109 
13043 re) 15.66287477 Rate — een aine 
10 15.66745489 15.6681 544 -0.00069951 
30 15.67900804 15.6788227 0.00018534 
50 15.68914747 15.6896365 -0.00048903 
100 15. 71437096 15.7173072 -0.00293624 
i 7840 0 tpcb6WNOOs?  ~' waaztacees. Sek cece 
10 14.56459216 14. 5646204 -0.00002824 
29 14. 56491839 14.5650397 -0.00012131 
52 14.56555106 14.5655565 -0.00000551 
103 14. 56649157 14. 5668036 -0.00031203 
wo 
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TABLE 4.8 


% MEAN MOTION VALUES FOR MODEL N2 
SATELLITE tt ACTUAL n CALCULATED nt” ERROR 
days rev/day rev/day rev/day 
15363 0 pT se > a ea 
10 1462055640 146183525 0.00220390 
30 14, 62850624 146310084 ~0.00250216 
50 14, 63611823 14, 6442566 -0.00813837 
100 1465988312 146799675 -0.01979188 
W760” O° SHSwrOIOHG: “uecceulce wens ane 
10 1447959223 14.4745202 0.00553513 
31 14,,48005842 144745202 0.00553822 
52 14,.48061500 14.4749875 0.00562750 
102 1448162641 14.4751171 0.00550931 
1303. 2~CO*t*=<‘“«~* G0 Metheny a ace ee 
10 15.66745489 15.6785692 ~0.01111431 
30 15.67900804 15.6884331 -0.00942506 
50 15. 68914747 15.6984351 -0.00928763 
100 15.71437096 15. 7240441 -0.00967314 
; 7840 0 UipSOULOOR2 eeczcecsue 7 ecto 
ie 10 14, 56459216 14.5602155 0.00437666 
29 14, 56491839 14, 5606034 0.00431499 
52 1456555106 14, 5610904 0. 00446066 
103 14, 56649157 14, 5622384 0.00425317 
oes 
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TABLE 4.9 


ACTUAL n 


14,.61592219 
1462055640 
14,.62850624 
1463611823 
1465988312 


1447940319 
14.47959223 
14,.48005842 
1448061500 
14.48162641 


15.66287477 
15.66745489 
15.67900804 
15.68914747 
15.71437096 


14. 56440052 
14.56459216 
14. 56491839 
14, 56555106 
14. 56649157 


MEAN MOTION VALUES FOR MODEL N3 


CALCULATED n 
rev/day 


14.6190976 
14.6320271 
14.6449565 
14..6772802 


144746000 
144750703 
144755406 
14..4766604 


15.6779276 
15.6879566 
15.6979856 
15.7230580 


14. 5606685 
14.5610645 
14, 5615438 
14. 5620670 


0.00145880 
-0.00352086 
-0.00883827 
-0.01739708 


0.00492230 
0.00498812 
0.00050744 
0.00496601 


-0.01047271 
-0.00894796 
-0.00883813 
-0.00868704 


0.00392366 
0.00385389 
0.00400726 
0.00442457 
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TABLE 4.10 


14.61592219 
1462055640 
1462850624 
14. 63611823 
14. 65988312 


14.47940319 
14.47959223 
14.48005842 
14.48061500 
14.48162641 


15.66287477 
15.66745489 
15.67900804 
15.68914747 
15.71437096 


1456440052 
14, 56459216 


14. 56491839 
1456555106 
14, 56649157 


MEAN MOTION VALUES FOR MODEL N4 


14.5997121 
146182269 
14.6361425 
146812678 


14.4470907 
14.4662946 
14.4851002 
14.5303213 


15.6907502 
15.7087211 
15.7269894 
15.7722919 


14. 5358578 
145529297 
14.5740243 
14.6201866 


0.02084430 
0.01027934 
. 00002427 
05138468 


03004711 
.02971306 
-0.03784193 
+ 05792094 


0.02873436 
0.01198869 
-0.00847324 
-0.05369503 
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TABLE 4.11 


14.61592219 
1462055640 
14.,62850624 
14, 63611823 
14..65988312 


14.47940319 
14.47959223 
14..48005842 
1448061500 
14.,48162641 


15.66287477 
15.66745489 
15.67900804 
15.68914747 
15.71437096 


14. 56440052 
14, 56459216 
14, 56491839 
14, 56555106 
14, 56649157 


MEAN MOTION VALUES FOR MODEL N5 


www we ew ew ee ww ee ow ww ww ew we = Ow ww ww on www a a aw en nn ee enw ee em mem em ewe ee 


CALCULATED n 
rev/day 


146225296 
14.635 7444 
14,6489592 
146819962 


wen eee ee 


14.4796342 
144801193 
14,.4806044, 
14.4817594 


15.6684158 
15.6794978 
15.6905798 
15.7182848 


14. 5646807 
14,.5652131 
14. 5658576 
14. 5672866 


ERROR 


-0.00198320 
~0.00723816 
~0.01284097 
-0.02211308 


-0.00004197 
-0.00006088 

0.00001060 
-0.00013299 


-0.00096091 
-0.00048976 
-0.00143233 
-0.00391384 


-0.00021484 
-0.00029471 
~0.00030654 
~0.00079503 


Since the criteria for quality of a model was set at 0.01595 


rev/day at 100 days, the absolute value of the error at 100 days 


(as near to 100 days as the data allowed) was averaged over the four 


satellites for each model. 


model quality and compare the models to each other. 


This average error was used to measure the 


Table 4.12 compares 


the 100 day error values of all satellites and all models. 
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0.41316866 
0.02909878 
0.0197918€ 
0.01739705 
0.05138468 
0.02211308 


TABLE 4.12 


0.00271859 
0.00024109 
0.00550931 
0.00496601 
0.04869489 
0.00013299 


0.09020714 
0.00293624 
0.00967314 
0.00868704 
0.05792094 
0.00391384 


100 DAY ERROR ABSOLUTE VALUES FOR ALL SATELLITES AND MODELS 


0.01188053 
0.00031203 
0.00425317 
0.00442457 
0.05369503 
0.00079503 


Table 4.13 shows the average of the 100 day error values for each 


model in order from smallest to largest. 


TABLE 4.13 


AVERAGE 100 DAY ERROR VALUES FOR EACH MODEL 


0.0067387 
0.0081470 
0.0088687 
0.0098068 
0.0529239 
0.1295098 


Analysis of the data for each model indicates the following: 


1) Models nO and n4¥ have a 100 day error greater than 0.01595 


rev/day and are removed from consideration. 
2) Models ni, n2, n3, and n5 all have acceptable 100 day error 


values as shown in table 4.14. 
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TABLE 4.14 


AVERAGE 100 DAY ERROR VALUES FOR ACCEPTABLE MODELS 


wee ee ee ew ee oe wn we we we we ew ee nn em ow = = = nw wn a ww wn a ww a a ee ee ee ee ee 


MODEL ERROR rev/day 
n5 0.0067387 
ni 0.0081470 
n3 0.0088687 
n2 0.0098068 
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3) Model n5 is the best in terms of 100 day error. It also 
produces very accurate values of n at times other than 100 days. 

4) Models n2 and n3 are very consistant in the magnitude of 
the error for any satellite as shown in Tables 4.8 and 4.9. This is 
shown in figures 4.1 and 4.2 which are plots of actual and calculated 
values on n. 

In studying the data it was found that the greatest error, for 
example for satellite 15363, occured when n was unstable. To eliminate 
this, N/2 for satellite 15363 was averaged over a 30 day span using 10 
day intervals. This decreased the 100 day error for that satellite 
by a factor of seven in the n5 model. The 30 day time span is indicated 
by the about 27 day length of the solar cycle. Since this first test 
improved the results a general test was done using equations 2.1 and 
2.2 for nand n. The times were chosen to be at the beginning, middle, 
and end of the 30 day period. Averaging over this time seemed to 
eliminate some of the randomness in the flux dependent values. The 


results are shown in tables 4.15 through 4.19. 
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Figure 4.1 Actual and Calculated values of n for Satellite 7840 using 


oe 


Unaveraged n and n 


TIME 
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Figure 4,2 Actual and Calculated values of n for Satellite 7840 using 


Averaged n and i 


wo 


0.00044071 
0.00001892 
0.00049658 
0.00001943 


TABLE 4.15 


AVERAGE VALUES OF n AND ni FOR EACH SATELLITE 


0.00000295 
0.00000029 
0.00000552 
0.00000022 


0.00023550 
0.00001156 
0.00026520 
0.00001057 


TABLE 4.16 


MEAN MOTION VALUES FOR MODEL Ni 


¥ 


COO GO RPE Ie PCR CR OPEL 


1462860624 
14, 63270865 
1464873275 
14.671 74206 


14.48005842 
14.48028419 
14.48126626 
14..48210420 


15.67900804 
15.68473145 
15.70452036 
15.73148592 


14. 56491839 
14.,56516626 
14, 56622862 
14.56703165 


1463320080 
14, 65203730 
14.67571370 


1448023990 
14.48109410 
14.,48212170 


15.68482970 
15.70557610 
15.73241580 


14. 56514620 


1456599960 
1456711360 


-0.00049215 
-0.00330445 
-0.00397164 


0.00004429 
.00017316 
.00001750 


- 00009825 
-00105574 
- 00092988 


.00002006 
- 00022902 
-00008195 
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yO TABLE 4.17 
eo oe MEAN MOTION VALUES FOR MODEL N2 
= SATELLITE ¢t ACTUAL n CALCULATED n—==~=*«CERROR 
By days rev/day rev/day rev/day 
* wwe ee Be wm ewe ew Dee ew ew OM HO em mw ww Be mM wee wm Om ew MO OM Mw eww wm ew Bem ew ewe we ew ew ew ee eee em ew we oe ewe eee eee 
15363 0 14.62860624 = -wnnnnee---- == +--+ -- 
10 14..63270865 1462934390 0.00336475 
50 14, 64873275 14. 64675540 0.00197735 
100 14.67174206 14. 66864400 0.00309806 
flee. 9° = Ge. TMBOOSBGS: " cnceeaswecs  ~ atesaaeewes 
9 14.48028419 14.47467020 0.00561399 
| 51 1448126626 1447545890 0.00580736 
iM 101 14,.48210420 1447641010 0.00569410 
Be 13043 0 1567000004. caeetcceee 7 en ewsensene 
* 11 15.68473145 15.69542640 -0.01069495 
50 15.70452036 15.71460470 ~0.01008434 
100 15.73148592 15.73942180 -0.00793588 
% 7840 0 SICENEIRSG.  -Cokececeeme token 
.. 11 14.,56516626 14. 56074740 0.00441886 
i 52 14. 56622862 1456153640 0.00469222 
1 & 105 14.56703165 14. 56256640 0.00446525 


TABLE 4.18 


14.62860624 
14.,63270865 
14..64873275 
14.67174206 


14.48005842 
14.48028419 
14.48126626 
14.48210420 


15.67900804 
15.68473145 
15.70452036 
15.73148592 


14.56491839 
14.56516626 
14. 56622862 
14.56703165 


MEAN MOTION VALUES FOR MODEL N3 


CALCULATED n 
rev/day 


14.62978980 
14.64744050 
14.66950400 


14.47521070 
14..47600930 
14.47696010 


15.69473160 
15.71415990 
1573906790 


14.56119420 


14.56199330 
14. 56302630 


rev /day 


- 00291885 
.00129225 
. 00223806 


00507349 
00525696 
00514410 


.01000015 
- 00963954 
. 00758198 


00397206 
. 00423532 
0.00400535 
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TABLE 4.19 


ie 
ad MEAN MOTION VALUES FOR MODEL N5 
SATELLITE Ne ACTUAL n CALQULATED r=——SCERROR 
days rev/day rev/day rev/day ' 
15363 0 Th 6pB6062h.. <Geesendee- 8 9 Sue eco 
10 14,63270865 1463221620 0.00049245 
50 14..64873275 14.,65205620 -0.00332345 
100 14.,67174206 14,67560620 ~0.00386414 
176 Ota NB0Seua  eceecseees | ete seesaee 
9 14.48028419 14.48026650 0.00001769 
51 1448126626 14.48123750 0.00002876 
101 1445210420 14.48239350 -0.00028930 
13043 > 0 15.67900804 eset = eee ea ‘ 
11 15.68473145 15. 68484240 -0.00011095 : 
50 15.70452036 15.70552800 -0.00100764 
100 15.73148592 15.73204800 -0.00056208 
7840 O° GAleololego. “\aicceetoens ~  * secueseetee 
11 14, 56516626 14.56515090 0.00001 536 5 
52 14, 56622862 14.56601770 0.00021092 
te 105 1456703165 14.,56713810 -0.00010645 
TABLE 4.20 


100 DAY ERROR ABSOLUTE VALUES FOR ALL SATELLITES AND MODELS 


MODEL/SATELLITE 15363 14476 13043 7840 
ne 0.00049215 0.0000175 0.00092988  0.00008195 
n2: 0.00309806 0.0056944 0.00793588 0.00446525 
n3: 0 .00223806 0.0051441 0.00758198 0.00400535 
n5: 0.00386414 0.0002893 0.00056208 0.00010645 
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TABLE 4,21 


AVERAGE 100 DAY ERROR VALUES FOR AVERAGED n AND n 


wee ee ee ee eee ee ee enw wen a a a ee ew ee ee a ee ee 


MODEL ERROR rev/day 
n5 0.0012055 
n3 0.0047424 
nl 0.0050009 
n2 0.0052980 


These values are about a factor of two (at least) improvement 
over the values produced by unaveraged n values. The NORAD model (n5) 
shows the most improvement and is still the most accurate. Its greatest 
problem is the instability of the error for a single satellite over 
time as shown in Table 4.19. In contrast, model n3, which has a small 
average error, shows great consistancy in the magnitude of the error 
as shown in Table 4.18. This means that if a value of n is known at 
some time around the start time, the size of the error to be expected 
from the model can be determined for that satellite. This could be of 
great help in locating lost objects. For example, if an elset is 
available for the time at which the satellite went lost, or start time, 
and two to three elsets are available for ten to thirty days prior to 
the start time, n and n may be calculated using the elset prior to the 
start time. Then a calculated value of n may be found and compared 
to the actual value of n at the start time to find a limit on the 
magnitude of the error. Next, at some future time, the value of n 
and a plus or minus range can be found. This will give the analyst 
general bounds within which to look for the satellite. For example 


satellite 519 has the following values of n at the times shown. 


4.22 


TABLE 4.22 


SAMPLE VALUES OF MEAN MOTION 


15.0054437 

15.0059095 

15.0062630 

15.0063419 

15.0091569 
Using equations 2.1 and 2.2 for times 34312, 84324, 84331 gives: 

n = 0.00003857 fm = 0,00000911 

Predicting from 84331 to 84338 using model n3 gives: n = 15.00824 
with an error of -0.0018981. Therefore the approximate error in the 
future is plus/minus 0.002. The actual error at 85072 is -0.0026. 
The plus/minus error gives the analyst a range within which to look for 
the satellite. 


An attempt was made to calculate n from a simple statistical 


model: 


yF (4.1) 


n(t) = Bo + Byn(t <1) 45 


+ 
at Bon, + 2B 
This model would be used in place of equation 2.1 to find n. BMDP was 
unable to produce an accurate result. The residuals from the model 


were of the same order of magnitude as the value being calculated. 


This results in no gain over conventional methods. 


V. CONCLUSIONS AND RECOMMENDATIONS 


The n5 model has been shown to be the most accurate when used 


with 30 day average values of N/2. The n3 model is a close second 


to the n5 model. Since the n3 model is stable, if some data for n 
near the the start time is available, limits on the error inn at 
some future time can be set. Because of this, the n3 model is better 
for the satellite orbit types studied here, that is, low eccentricity 
orbits with 14 n 16. 

Further study is recommended in four main areas. The first area 
for further study is the possible use of N/2 to find a value for n 
and the use of these values in model n3. This would determine if the 
greater accuracy available in the n5 model can be combined with the 
consistancy of the n3 model to give more accurate in-track areas in 
which to search for lost objects. 

A second item for further study is the extension of the model to 
times greater than 100 days for prediction. 

The third subject for study is the inclusion of mean motions 
outside the range included in this thesis and for more eccentric 
satellite orbits to be included. 

The final study should determine if other orbital elements need 


to be included in the model. 
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n(t) 


0.0000401 
0.9000286 
0.90000505 
0, 0000558 
0.0000668 
0, 9000849 
0O.9000629 
0.0000458 
oO,0000525 
0.0000615 
oO, 0000690 
0, 0000854 
o.900082S5 
9.9001192 
0.9000792 
O.9000563 
0.0000896 
0, 0000656 
0. QO000150 
O. 9000208 
0. 9000334 


0, 90002038 
O.9000149 
0.0000222 
O,0000101 
0O.,0000105 
O.90000146 
O.0000141 
0.0000207 
9.90000211 
0, 00002253 
0,.0000181 
Oo.0000134 
0.9000105 


APPENDIX A 


DATA FILE USED FOR BMDP 


n(t - dt) 


Oo, O000000 
O.0000401 
Oo. 0000286 
0, OOOOSOS 
Oo, OOO0SSB 
0. 0000668 
O.0000849 
0.0000629 
o.0000458 
0.0000525 
0.0000615 
0.0000690 
0. 0000854 
0, 0000825 
O.9001192 
0. 0000324 
O.0000792 
0. 0000563 
0. 0000896 
Oo, 0000000 
O,0000150 
0.0000208 
0, O000334 
0,.0000203 
0.0000149 
O, OOOO222 
O.Q00010L1 
O.0000105 
Oo. Q0O00000 
0O,0000141 
QO, 0000207 
O.0000211 
0,0000223 
0o.0000181 
0O.0000134 


n(t) 


14. 9814358 
14.9814758 
14. 9815617 
14. 9817133 
14. 9822159 
14. 9825630 
14. 9827328 
14.9829216 
14.9829674 
14. 9830725 
14. 9833183 
14. 9835253 
14..9826960 
14. 9838610 
14. 9839802 
14. 9840450 
14.9841242 
14. 9842367 
14. 9844160 
14. 8682432 
14. 8683929 
14.8685179 
14.86855135 
14. 8686934 
14.8688126 
14. 8689012 
14.8689518 
14. 86896235 
14.6311455 
14.6312017 
14,.6312627 
14, 6313906 
14.6316137 
14.6316862 
14.62317263 


A.1 


n(f) 


14.9813442 
14.9812442 
14.9813442 
14.9813442 
14.9613442 
14.9813442 
14.9813442 
14. 9812442 
14.9813442 
14.9813442 
14.9813442 
14.9813442 
14.9813442 
14.9813442 
14.9813442 
14.9813442 
14.9813442 
14. 9813442 
14.9813442 
14. 8681669 
14.8681669 
14.8681669 
14.8681669 
14.8681669 
14.8681669 
14.8681669 
14.8681669 
14.8681669 
14.6310959 
14.6310959 
14.6310959 
14. 6310959 
14.6310959 
14.6310959 
14.62310959 
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0.0000160 
0.00001 41 
O,0000111 
0, 0000312 
0, 0000238 
QO, OOOO300 
0. 0000243 
0,0000606 
0O.9000391 
0. 9000982 
Oo. 00004358 
0, 0000367 
0, 0000591 
0. 9000248 
0,0000264 
o,9000200 
0.90000267 
O,0000341 
0. 0000277 
O.0000404 
0,0000715 
0,0000467 
0. 9005770 
0,.0007029 
ao.90012178 
0.0012465 
0.0005617 
0,00073572 
Q.,0010633 
0.9010376 
0.0016823 
9.0014992 
0,0017462 
0.9015879 
0,.0014381 
0.0017490 
0O.0017366 
0,.90016155 
.0017204 


0.0000105 
0.90000160 
0.00001 41 


Q0,0000312 
0.0000238 
0, 0000300 
0.0000243 
0.0000606 
0,0000391 
Q, 0000982 
0o,0000458 
©. 0000367 
0,0000591 
0.0000248 
O.0000264 
0.0000200 
0.0000267 
O,0000241 
0,0000277 
0,0000404 
0,0000715 
O.0000467 
OoO.9005770 
0,0007029 
0,.0012178 
0.90012465 
0, 00035617 
0.0007572 
0.0010633 
0,9010376 
0.00168235 
0,.9014992 
0.9017462 
0.0015879 
0O.0014381 
0,0017490 
0,.0017366 
O.,90016155 


14.6317577 
14.6318216 
14.6318922 
15.0053215 
15.0054464 
15.0054941 
15.0055542 
15.0056515 
15.0057726 
15.0058117 
15. 0059891 
15. 0060806 
15.0061541 
15.0062132 
15.0062628 
15.0063419 
15. 0063820 
15.0064354 
15.0065718 
15.0966824 
15. 0068025 
15.0068750 
15.3779602 
15.3789206 
15.3796234 
15. 38084153 
15.3831081 
15. 3840485 
15. 3855299 
15.3875751 
15.3899431 
15.3922529 
15.3937531 
15.3957367 
15. 3973246 
15.3987627 
15.4005117 
15. 4022484 
15.4046593 


A.2 


14.63510959 
14.6210959 
14.6210959 
15.0049667 
15.0049667 
15.0049667 
15.0049667 
15. 0049667 
15.0049667 
15.0049667 
15. 0049647 
15.0049667 
15.0049667 
15.0049667 
15.0049667 
15.0049667 
15.0049667 
15.0049667 
15.0049667 
15.0049667 
15.90049667 
15.0049667 
15.3755236 
15. 3755236 
15. 3755236 
15. 3755236 
15. 3755236 
15. 3755236 
15. 3755236 
15. 3755236 
15. 3755236 
15. 3755236 
15. 2755236 
15. 3755236 
15. 3755236 
15. 3755256 
15. 2755236 
15.3755236 
15. 3755235 
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18 1 O.0013862 0.0017204 15.4063797 15.37552356 i -1 74 
20 = 0.0026684 0,0013862 15.4091520 15.3755236 -27 1 74 
= 1 O.0032047 0,0026684 15.4118204 15. 3755236 1-3 75 
2 3 0.9019274 9,0032047 15.4214344 15.3755236 6 1 73 
25 1 0.0006914 0.0019274 15.4243202 15.3755236 o © 74 
2 1 0.0007496 0,0006914 15.4255304 15. 3755226 Oo oO B80 
2 1 O.0011072 0.0007496 15.4262800 15.3755236 1 © 80 
28 1 0.0018005 0.0011072 15.4283867 15.3755236 2 2 80 
2 1 9,0020504 0.0018005 15.4307737 15.3755236 2 2 81 
30 1 0.0022745 0,0020504 15.4328241 15.3755236 Oo 2 8F 
31 1 0.0017424 0.0022745 15.4350986 15.3755236 1 © 85 
32 1 0.0013084 0.0017424 15.4372549 15.3755236 -3 -3 85 
33 1 0.0025711 0.0012084 15.4393463 15.3755236 ~-4 -4 86 
34 1 0.0018959 0.0025711 15.4419174 15.3755236 Oo -4 83 
35 1 9.0021210 0.0018959 15.4438133 15.3755236 Oo © 79 
36 1 O,.0021133 0.0021210 15.4476538 15.3755236 0 0 79 
38 2 0.0016947 0.0021133 15.4548283 15.3755236 ~-2 -2 79 
ae 40 2 0.0012426 0.0016947 15.4597807 15.3755236 -1 -1 78 
l¢ 41 i 9.0012321 0.0012426 15.4616613 15.3755236 Oo 0 76 
42 1 0.0024109 0,0012321 15.4628955 15.3755236 1 © 75 
3 1 ©.0017233 0.0024109 15.4653044 15.3755236 1 1 75 
44 1 0.0010376 0.0017233 15.4681396 15.3755236 1 1 76 
45 1 0.0005436 0,0010376 15.4779492 15.3755256 3.3 77 
46 1 0.0025892 0.0005436 15.4714928 15.2755256 -1 3 78 
47 1 0,0013638 0.0025892 15.4756136 15.3755236 ~-1 -1 81 
48 1 0.0022907 0.0013638 15.4769773 15.3755236 ~-2 -1 80 
49 1 0.0012608 ©.0012608 15.4815140 15.3755256 Oo Oo 79 
so 1 0.0035906 0.0012608 15.4827747 15.3755236 Oo OO 77 
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APPENDIX 8 


NORAD ELEMENTS FOR SATELLITES 15363, 14476, 13043, 7840 
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